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Abstract: We prove a change of variable formula for the 2D fractional Brownian motion of index 
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1 Introduction and main result 

In [I], Coutin and Qian have shown that the rough paths theory of Lyons [13] can be applied to the 
2D fractional Brownian motion B = (B^\ B^) under the condition that its Hurst parameter H 
(supposed to be the same for the two components) is strictly bigger than 1/4. Since this seminal 
work, several authors have recovered this fact by using different routes (see e.g. Feyel and de La 
Pradelle [7], Friz and Victoir [8] or Unterberger [19] to cite but a few). On the other hand, it is 
still an open problem to bypass this restriction on H . 

Rough paths theory is purely deterministic in essence. Actually, its random aspect comes 
only when it is applied to a single path of a given stochastic process (like a Brownian motion, 
a fractional Brownian motion, etc.). In particular, it does not allow to produce a new alea. As 
such, the second point of Theorem II. 21 just below shows, in a sense, that it seems difficult to reach 
the case H = 1/4 by using exclusively the tools of rough paths theory. 

Before stating our main result, we need some preliminaries. Let W be a standard (ID) 
Brownian motion, independent of B. We assume that B and W are defined on the same probability 
space (Q,^,P) with & = cr{B} V <t{VF}. Let (X n ) be a sequence of cr{-B}-measurable random 

variables, and let X be a ^"-measurable random variable. In the sequel, we will write X n s *^? x 
if (Z,X n ) (Z,X) for all bounded and cr{I?}-measurable random variable Z. In particular, we 
see that the stable convergence imply the convergence in law. Moreover, it is easily checked that 
the convergence in probability implies the stable convergence. We refer to [TI] for an exhaustive 
study of this notion. 

Now, let us introduce the following object: 
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Definition 1.1 Let f : R 2 — ► R be a continuously differentiable function, and fix a time t > 0. 
Provided it exists, we define J* * V f(B s ) ■ dB s to be the limit in probability, as n — > oo, of 



Lntj-l dl(p>m o(2) \ , S£/ R (l) o(2) 

dx \ D k/n' D kln> dx ^(fc+D/n' -°fc/n _ 

"(fc+l)/n ~~ ^k/n) 



T / . \ Ox\- k/n' k/n> ^_ dx\ ±J (k+l)/n' k/n' , p (l) R (l) 

'"W :_ 2^ 9 l-°(fe+l)/n ~ , 

fc=0 



Lntj-l o(2) n , df( B m o(2) x 

d y y D k/n' n k/n> ~T ch/^fc/n'^fc+iyJ / (2) R (2) s nn 

+ 9 l-°(ifc+l)/n -°fc/n^ 



fc=0 



If I n (t) defined by fH2P ^oes no£ converge in probability but converges stably, we denote the limit 
by f*Vf(B s )-d*B s . 

Our main result is as follows: 

Theorem 1.2 Let f : M? — > R be a function belonging to 'rf 8 and verifying (Hg), see 113.15)) below. 
Let also B = (B^ 1 ' , B^) denote a 2D fractional Brownian motion of Hurst index H £ (0, 1), and 
t > be a fixed time. 

1. If H > 1/4 then Vf(B s ) ■ dB s is well-defined, and we have 

f(B t ) = /(0)+ /' t Vf(B 8 )-dB s . (1.2) 
Jo 

2. If H = 1/4 then only £ Vf(B s ) ■ d*B s is well-defined, and we have 

f(B t ) L =* /(0) + £ Vf(B s ) ■ d*B s + ^j=f ^( B s) d Ws. (1-3) 

Here, cr^ is the universal constant defined below by 111.5]) . and f td2 f/dxdy(B s )dW s denotes 
a classical Wiener integral with respect to the independent Brownian motion W . 

3. If H < 1/4 then the integral J Q * B s -d*B s does not exist. Therefore, it is not possible to write 
a change of variable formula for B^ B^ using the integral defined in Definition \l.ll 

Remark 1.3 1. Due to the definition of the stable convergence, we can freely move each 
component in (|1.3p from the right hand side to the left (or from the left hand side to the 
right ) . 

2. Whenever (3 denotes a one-dimensional fractional Brownian motion with Hurst index in 
(0, 1/2), it is easily checked, for any fixed t > 0, that X^L=d 1 Pk/n{P(k+i)/n ~ Pk/n) does 
not converge in law. (indeed, on one hand, we have 

\nt\-l \nt\~l \nt\-l 

P\nt\/t= ^2 (P(k+l)/n ~ Pk/n) = 2 ^2 (P(k+l)/n ~ Pk/n) + ^2 {P(k+l)/n ~ Pk/n) 

k=0 k=0 k=0 

and, on the other hand, it is well-known (see e.g. [T2]) that 

n 2H-l (%+!)/„-&/„) 
k=0 
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These two facts imply immediately that 

\nt]-l / \nt\-l 



E Pk/n{P(k+l)/n ~ Pk/n) ~ 9 ( Pfntj/t ~ E {P(k+l)/n ~ Pk/n)' 



k=0 ~ \ k=0 



does not converge in law). On the other hand, whenever H > 1/6, the quantity 

[nij-l 

E ^{fiPk/n) + /(/3(fc+l)/n)) (/3(Jfc+l)/n - Afc/n) 
fc=0 

converges in L 2 for any regular enough function / : R — ► M, see [9] and [3]. This last fact 
roughly explains why there is a "symmetric" part in the Riemann sum (jl.lj) . 

3. We stress that it is still an open problem to know if each individual integral J* * ^(B s )d^Bs 
and Jq {gt(B s )d^ Bg 2 ^ could be defined separately. Indeed, in the first two points of Theo- 
rem 03 we "only" prove that their sum, that is f*Vf(B s ) ■ d^B s , is well-defined. 

4. Let us give a quicker proof of (|1.3|) in the particular case where f(x, y) = xy. Let f3 be 
a one-dimensional fractional Brownian motion of index 1/4. The classical Breuer-Major's 
theorem [T] yields: 

[n-J-l [n-J-1 

4= E (^(% + D/n-A/n) 2 -i) L = E (03H-i-A) a -i)H^w- d- 4 ) 

V fc=0 V fc=0 

Here, the convergence is stable and holds in the Skorohod space S> of cadlag functions 
on [0,oo). Moreover, W still denotes a standard Brownian motion independent of (3 (the 
independence is a consequence of the central limit theorem for multiple stochastic integrals 
proved in [18]) and the constant 0*1/4 is given by 



01/4 



\ J] + + - 2v4^) 2 < 00. (1.5) 
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Now, let /3 be another fractional Brownian motion of index 1/4, independent of [3. Prom 
lfL~4"|) . we get 

\nt\-l \nt\-l 



—7= E {v^(P(k+l)/n ~ Pk/n) 2 ~ l)) —r= E (V^(^(fc+l)/n ~~ A/n) 2 ~ l) 
V n t_n V n ,„_ n 

01/4 (W,WQ 



fc=0 v fc=0 

stably 



for (W, W) a 2D standard Brownian motion, independent of the 2D fractional Brownian 
motion (f3,{3). In particular, by difference, we have 

[n-J-l 

\ E ((%+l)/n-/Vn) 2 "(%^ 

fc=0 ^ 
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Now, set £?W = (j3 + (5)/V2 and B^ = (/3-/3)/^2. It is easily checked that flt 1 ) and B^ 
are two independent fractional Brownian motions of index 1/4. Moreover, we can rewrite 
the previous convergence as 

e'c. - - <t) h 7f * (»■«) 



fc=0 



with B^ l \ B^ and independent. On the other hand, for any a,b,c,d £ R: 
bd — ac = a(d — c) + c(6 — a) + (6 — a)(d — c). 

Choosing a = flg^ 6 = B^ +1)/n , c = Bfj n and d = £ ( ( ^ 1)/n , and suming for k over 
0, . . . , |_77.t J — 1, we obtain 

[ntj-l 

ijf 1 ) Bp) = Y fl« (Bfl 1V - Bw ) + B« (B« - Bf» ) 

\nt\/n \nt\/n / j k/n\ [k+l)/n k/nl k/n\ (k+l)/n k/nJ 

k=0 
\nt\-l 

+ E ( B Sil)/n " O «i)/n - O • 
fc=0 



Hence, passing to the limit using (|1.6|) . we get the desired conclusion in (|1.3p . in the par- 
ticular case where f(x,y) = xy. Note that the second term in the right-hand side of (11.71) 
is the discrete analogue of the 2-covariation introduced by Errami and Russo in [6]. 

5. We could prove (|1.3p at a functional level (note that it has precisely been done for f(x, y) = 
xy in the proof just below). But, in order to keep the length of this paper within limits, we 
defer to future analysis this rather technical investigation. 

6. In the very recent work [16], Reveillac and I proved the following result (see also Burdzy 
and Swanson [2j for similar results in the case where (5 is replaced by the solution of the 
stochastic heat equation driven by a space/time white noise). If j3 denotes a one-dimensional 
fractional Brownian motion of index 1/4 and if g : R — > M. is regular enough, then 

- T Y,9(Pk/n){MP(k + l)/n-Pk/n) 2 -l) S ^l l / fl" \Ps)ds + <7 1/4 / g(J3 s )dW, (1.8) 



k=0 

for W a standard Brownian motion independent of f3. Compare with Proposition 13. 3l below. 
In particular, by choosing g identically one in (|1.8p . it agrees with (|1.4D . 

7. The fractional Brownian motion of index 1/4 has a remarkable physical interpretation in 
terms of particle systems. Indeed, if one consider an infinite number of particles, initially 
placed on the real line according to a Poisson distribution, performing independent Brownian 
motions and undergoing "elastic" collisions, then the trajectory of a fixed particle (after 
rescaling) converges to a fractional Brownian motion of index 1/4. See Harris |10| for 
heuristic arguments, and Diirr, Goldstein and Lebowitz [5] for precise results. 

Now, the rest of the note is entirely devoted to the proof of Theorem 11.21 The Section 2 
contains some preliminaries and fix the notation. Some technical results are postponed in Section 
3. Finally, the proof of Theorem 11.21 is done in Section 4. 
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2 Preliminaries and notation 



We shall now provide a short description of the tools of Malliavin calculus that will be needed in 
the following sections. The reader is referred to the monographs [14] and [17] for any unexplained 
notion or result. 

LetS = (sfUf )te[o,T] be a 2D fractional Brownian motion with Hurst parameter belonging 
to (0, 1/2). We denote by TL the Hilbert space defined as the closure of the set of step ]R 2 -valued 
functions on [0, T], with respect to the scalar product induced by 

((^■[o,*!]! 1 ^])) { 1 [o,s 1 ]A[o,s 2 ])) n = RH{h,si) + RH(t 2 ,s 2 ), Si,U G [0,T], i = 1,2, 

where R H (t, s) = \ (t 2H + s 2H - \t - s\ 2H ). The mapping (l [0 , tl ], l[ ,t 2 ]) ^ + B^ can be 
extended to an isometry between TL and the Gaussian space associated with B. Also, f) will 
denote the Hilbert space defined as the closure of the set of step R- valued functions on [0, T], with 
respect to the scalar product induced by 

( ![<),*]) 1 [0,a]) 5 = R H (t,s), s,te [0,T]. 

The mapping l[ 0j t] •— ► 5* equals 1 or 2) can be extended to an isometry between 9) and the 
Gaussian space associated with B^\ 

Consider the set of all smooth cylindrical random variables, i.e. of the form 

F = f(B(<p 1 ),...,B(<p k )), (pi 6 Ti, i = l,...,k, (2.9) 

where / G < ^' 00 is bounded with bounded derivatives. The derivative operator D of a smooth 
cylindrical random variable of the above form is defined as the W-valued random variable 

DF = J2jf-. i B &i), • • • » B(<Ph))<Pi =■ {D BW F, D BW F) . 



dxi 
i=i 



In particular, we have 



D B (i)Bi = Sijl\ 0>t ] for i,j G {1,2}, and 8^ the Kronecker symbol. 



By iteration, one can define the mth derivative D m F (which is a symmetric element of L (Q, H® m )) 
for m ^ 2. As usual, for any m ^ 1, the space ID m ' 2 denotes the closure of the set of smooth 
random variables with respect to the norm || • \\ m ,2 defined by the relation 



\f\\ 2 ^ 2 = e\f\ 2 + J2e\\d 1 f 



2 

1=1 



The derivative D verifies the chain rule. Precisely, if ip : W 1 — > R belongs to c tf 1 with bounded 
derivatives and if Fi, i = 1, . . . ,n, are in D 1,2 , then ip(F±, . . . , F Tl ) G O 1,2 and 

D<p(Fi, . . . ,F n ) = ■ ■ ■ , Fn)DFi. 

i=i 
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The mth derivative D™^ (i equals 1 or 2) verifies the following Leibnitz rule: for any F,G G 
such that FG G D m ' 2 , we have 



W«^* ll ..,U = E(^«^). 1 ,..^W^„ lt ... |Ilm _ f .. <i€[0,n i = l,... >TO> (2.10) 

where the sum runs over any subset {si, . . . , s r } C {t\, . . . , £ m } and where we write {t±, . . . , t m } \ 

{si, , s r } — . {u^,..., U m — r } . 

The divergence operator 5 is the adjoint of the derivative operator. If a random variable 
u G L 2 (Q,TC) belongs to domJ, the domain of the divergence operator, then 5(u) is defined by the 
duality relationship 

E(F5(u)) =E(DF,u) n 

for every F G ED 1,2 . 

For every q ^ 1, let Ti. q be the qth Wiener chaos of B, that is, the closed linear subspace of 
L 2 (Cl,A,P) generated by the random variables {H q (B (h)) ,h G Tt, \\h\\y, = 1}, where H q is the 
qth Hermite polynomial given by H q (x) = (— l) q e x ^ 2 -^{&~ x //2 )- The mapping 

I q (h^) = H q {B{h)) (2.11) 

provides a linear isometry between the symmetric tensor product TL &q (equipped with the modified 
norm ~^\\ ' Wn^i) an d ^<?- The following duality formula holds 

E(FI q (f)) = E((D q F,f) Hm ), (2.12) 

for any / G 7i Gq and F G ED 9 ' 2 . In particular, we have 

E(FiP( 9 ))=E(( D i m F, 9 )^), 4 = 1,2, (2.13) 

for any g G f) 09 and F G ID 9 ' 2 , where, for simplicity, we write Iq\g) whenever the corresponding 
gth multiple integral is only with respect to SW. 

Finally, we mention the following particular case (actually, the only one we will need in the 
sequel) of the classical multiplication formula: if /, g G f), q ^ 1 and i G {1,2}, then 

j(*)(/^)I«(^)=^r!^ (2.14) 

r=0 

3 Some technical results 

In this section, we collect some crucial results for the proof of (II. 3|) . the only case which is difficult. 
Here and in the rest of the paper, we set 

AB k/n := B (k+l)/n ~ B k)n> S k/n ■= l[fc/n,(fc+l)/n] and £ k/n l[0,fc/«]' 

for any % G {1, 2} and k G {0, . . . , n — 1}. 

In the sequel, for g : M? — > E belonging to < ^" ? , we will need assumption of the type: 



(H q ) sup E 

se[o,i] 



dx a dy' 



<bP, bp: 



< oo for all p ^ 1 and all integers a, 6 ^ s.t. a + b ^ g. 
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We begin by the following technical lemma: 
Lemma 3.1 Let (3 be a ID fractional Brownian motion of Hurst index 1/4. We have 



(i) \E(Pr(j3t - S )) | < y/\t -s\ for any < r, s, t < 1, 

n-l 

I At., A, , \ 

n, 



^ X] \( £ l/n^k/n 
k,l=0 

n-l 



fc,«=0 

n-l 



fc=0 
n-l 



( e fc/n> ^ 



k/n/ 



+ 



2^ 



( e fc/n; ^fc/n 

fc=0 

Proof of Lemma 13.11 

(i) We have 



2 _ J_ 
£ 4n 



0(n), 

O^ 1 "^ 2 ) for any r ^ 1, 
= O(l), 

>oo 

O(Vv^). 



- &)) = ^(v 7 * - v^) + J (vV 



vit" 

2 V v ' 2 V v 1 1 Vl 
Using the classical inequality \y/\b\ — \f\a\ \ < — of, the desired result follows. 



(ii) Observe that 



'\ £ l/ni$k/n) , 



1 



Vk+l-Vk- y/\k + l-l\ + y/\k 



Consequently, for any fixed I G {0, . . . ,n — 1}, we have 

n-l /z-i 

£|v^/nA/n)J ^ 2 + 2~7^ J^v 7 ^- ^-^-1 

fe=0 V \fc=0 

n-l \ 

+1+ Vfc - / + l - v 7 ^! 

fc=Z+l / 



n-l 



from which we deduce that sup N ( e z/n5 ^fc/n)^ = 0(1)- It follows that 



0<Z<n-l 



fc=0 



n-l 



n-l 



^ (£|/«>4/«> s sup ^2\(ei/ n ,6, 



k,l=0 



0<Z<n-l 



k/n/ 



0{n). 



(3.15) 



fc=0 
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(Hi) We have, by noting p(x) = +T| + yj\x — Tj — 2-y/jxj): 



n-l 



n-l 



n 



-r/2 



^l/ti-^Kn 1 -/^!^)!. 



fc,Z=0 fc,/=0 feGZ 

Since X^feez l/° r (^)l < oo if r ^ 1, the desired conclusion follows. 
(iv) is a consequence of the following identity combined with a telescopic sum argument: 



1 



2v^ 



2v^ 



fttj We have 



— f \/fcTl - VA; + 1 - v 7 ^ - 2 
4re V / 



Thus, the desired bound is immediately checked by combining a telescoping sum argument 
with the fact that 

1 



Vk + 1- Vk-2 



Vk + l + Vk 



< 2. 



□ 



Also the following lemma will be useful in the sequel: 

Lemma 3.2 Let a ^ and q ^ 2 be two positive integers, g : M 2 — > R be any function belonging 
to c € 2q and verifying (H2 q ) defined by $3.15\) . and B = (B^\ B^) be a 2D fractional Brownian 
motion of Hurst index 1/4. Set 

n-l 
k=0 

where H q denotes the qth Hermite polynomial defined by H q (x) = (— l) l? e a;2//2 ^- (e~ x2 ^ 2 ) . Then, 
the following bound is in order: 



E(\V n \ 2 )=0( 



l-q/2-a/2 



) as n — > oo. 



(3.16) 
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Proof of Lemma I3.2L We can write 



e (iki 2 ) e ^[^g..<M<.^(^S!J a ( a O q 



fc,Z=0 

n-l 



J27T11 



fc,Z=0 

« — n \ / i. n 



E * Q E E ("ST 



r=0 v 7 fc,Z=0 



r=0 V 7 a+b=2g-2r ' ' fc,/=0 V tf y 

x K)° (29 " 2r)! (^f/n* ^ 8 ^JT)^ fo/n, <«■ 

(3.17) 

Now, observe that, uniformly in k, I £ {0, . . . , n — 1}: 

€V ® «T ) „ _ t_ 0(n-("-)), see Lemma O (i) , 

= 0(n~ a / 2 ), use (H 2q ), 



■k/n^l/n^k/n & °l/n / m2q - 2r n Zoo 

E ( P-P?) ,Bf] )^(B« B#>)(A2$ r MM 

V dll a «/"• fc/W 'All I l/n'\ k/n/ \ l/n/ I 



and, also: 

n-l 

(Sk/n^i/nYsj = 0(n 1 ~ r / 2 ) for any fixed r ^ 1, see Lemma [37T] (iii) . 

fe,i=0 

Finally, the desired conclusion is obtained by plugging these three bounds into (|3.17|) . after having 
separated the cases r = and r = 1. □ 

The independent Brownian motion appearing in (|1.3j) comes from the following proposition. 

Proposition 3.3 Let (/?,/?) be a 2D fractional Brownian motion of Hurst index 1/4. Consider 
two functions g,g : ]R 2 — ► K belonging in 'rf 4 , and assume that they both verify (H4) defined by 
(MB- Then 

n—l ^ n— 1 

-p^9(A/n>A/n)(\^(%n) 2 ~ l) , ^ ^ (A/n, A/n) (\/^(A^ fc / n ) 2 - l) 

v n fc=o ^ n fc=0 

stably/ ^ /<0 ~, JTI7 1 /"^V* ~ ~. ~ 1 Z- 1 ^ 



n— +00 



9(Ps,0s)dW 9 + ±J yt(/3 s ,p s )ds,a 1/4 ^ g(0.,Ps)dW s + ±J ^(p s J s )ds\ 
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where (W,W) is a 2D standard Brownian motion independent of (f3,/3), and (T\u is defined by 

In the particular case where g(x,y) = g(x) and g(x,y) = g(y), the conclusion of the proposition 
follows directly from (jl.81) . In the general case, the proof only consists to extend literaly the proof 
of (11.81) contained in \W\. Details are left to the reader. 



4 Proof of Theorem 11.2 



We are now in position to prove our main result, that is Theorem 11.21 

Proof of the third point (case H < 1/4). Firstly, observe that (|1.4p is actually a particu- 
lar case of the following result, which is valid for any fractional Brownian f3 with Hurst index H 
belonging to (0,3/4): 

V k=0 

with W an independent Brownian motion and ajj > an (explicit) constant. By mimicking the 
proof contained in the fourth point of Remark 1.3, we get, here, for any H G (0,3/4), 



n 2ff-l/2 y 1 AB (1) AB (2) ^ (4 _ 18) 

k/n k/n n ^ oc ^ 



k=0 

But, see (jl.7p . the existence of f Q B s ■ d*B s would imply in particular that ^l=o 1 ^^k/n^^k/n 
converges in law as n — > oo, which is in contradiction with (|4.18p for H < 1/4. The proof of the 
third point is done. 

Proof of the second point (case H = 1/4). For the simplicity of the exposition, we assume 
from now that t = 1, the general case being of course similar up to cumbersome notation. For 
any a,b,c,d G R, by the classical Taylor formula, we can expand f(b,d) as (compare with (|1.7|) ): 

f{a, c) + dif(a, c){b - a) + 0*2/(0, c)(d - c) + -duf{a, c){b - a) 2 + ^3 22 /(a, c)(d - c) 2 

+ ldmf(a, c)(b - af + \d 222 f(a, c){d - cf + ^o\m/(o, c)(b - a) 4 + ±-d 2222 f(a, c)(d - c) 4 
b b 24 24 

+ 6\ 2 /(a, c)(b -a)(d-c) + ^d U2 f(a, c)(b - af(d - c) + \d 122 f{a, c)(b - a)(d - cf 

+ \d lll2 f(a, c){b - af(d - c) + \d ll22 f{a, c)(b - af(d - cf + \d l222 f(a, c){b - a){d - cf 
6 4 6 

(4.19) 

plus a remainder term. Here, as usual, the notation 6V..12...2/ (where the index 1 is repeated k 
times and the index 2 is repeated I times) means that / is differentiated k times w.r.t. the first 
component and I times w.r.t. the second one. By combining (|4. 19|) with the following identity, 
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available for any h : R — > R belonging to ^ : 

fc'(a)(6 - a) + ^"(a)(6 - a) 2 + jU'"(a)(6 - a) 3 + ^fc""(a)(6 - a) 4 
= — - - ^ — —(6 — a) — — h"'(a)(b — a) 3 — —h""{a){b — a) 4 + some remainder 
we get that f(b,d) can also be expanded as 

lo m /(a,c)(6-a) 3 -i 



/(a, c) + i (^/(a, c) + 9i/(6, c)) (6 - a) - ^d lu f(a, c)(b - a) 3 - ^dim/(a, c)(6 - a) 4 
1 ,„ , „ „, „ w . . 1 „ , 1 



+ ± (a 2 /(a, c) + 8 2 f(a, d)) (d-c)- -^d 222 f(a, c)(d - cf - ^8 2222 f(a, c)(d - c) 
+ di 2 f(a, c)(b - a)(d - c) + -<9ii 2 /(a, c){b - a) 2 (d - c) + -<9i 22 /(a, c)(6 - a)(d - c) 2 
+ ^112/(0, c)(6 - af(d - c) + j5 1122 /(a, c)(6 - a) 2 {d - c) 2 + \d 1222 f(a, c)(b -a)(d- c) 3 

D 4 D 

(4.20) 

plus a remainder term. 

By setting a = B^j n , b = B$ +1 y n , c = BfJ n and d = B^ +l)/n in (|4.20|), and by suming the 
obtained expression for k over 0, ... ,n — 1, we deduce that the conclusion in Theorem 11.21 is a 
consequence of the following convergences: 

n— 1 _ „i 

S§> := EWa^jKU'^- / W(B?\Bp>)* (4.21) 

j_ n J 



fc=0 
n-l 



E^iH/(</K)(A0 4 3 / W(B« (4-22) 

^ := EW*<)K) 3 ^- W(BM 2 >)* (4-23) 
fc=0 Jo 
n—l ,.\ 

■= £W(<,<)(A2^ 4 £-3/ d 2222 f(BP,BP)ds (4.24) 
fc=o 1/0 

fc=0 

+\j^d ll22 f{B^\B^)ds (4.25) 

n— 1 1 „i 

S<«> :_ E^/t^.B^JtA^J^BW^-I/ WWf.flf)* (4-26) 

n—l 1 „x 

^ := E B$,)AB« (Alf/j 2 £^~\ d 1122 f(BPM»)ds (4.27) 

J n w 
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n— 1 „i 

^ - EW(<,<)(A<;j 2 (A<J 2 n -£/ W(*<W)* (4.28) 

l n "'O 



fc=0 
n-l 



,(1)\3 A „(2) Prob 



fc=0 
n-l 



(4.29) 



(4.30) 



fc=0 



Note that the term corresponding to the remainder in (|4.20|) converges in probability to zero due 
to the fact that B has a finite quartic variation. 

Proof of UJT\} , (THity , (iMj) and W¥§ . By Lemma E21 with q = 3 and a = 0, and by using 
the basic fact that 



'fc/r 

we immediately see that (|4.23|) is a consequence of the following convergence: 
i n— i i pi 

4 £ *-/(<• £i - 5 / w<a? >, *?>)*. 

v fc=0 1/0 

So, let us prove (j4.32[) . We have, on one hand: 



E 



1 n— 1 

V fc=0 



(2) 
fc/n 



n-l 



k,l=0 
n-l 



= \ E E ( 5 222/(< n , Bf] n ) d 222 f(Bljl B\f n ) if fo /n ® fy B ) 

fc,«=0 
1 n— 1 

+ ~E E ( 5 222/(< n , Bg,) W(*#, 2J« )) (8 k/n , 5 l/n)S) 
k,l=0 
^ n— 1 

= ~ E ^ ( <9 22222/(^ 1 / ) n , £f/ } J d 2 22f(Bljl, B^J^j (s k/n , 5 k/n ) ft{s k/n , 5 l/n )sj 
k,l=0 
1 n— 1 

+ - E E (<W( S */L- B k/J 0*nHB% Bg)) <£*/„, «*/»>*<*,/„, 

fc,Z=0 
1 n— 1 

+ - E E { d ^f( B k/l B k/n) d ^f( B l/i B lfj) ^ Sk/n)*{ek/n, «!/»>« 
fc,«=0 

1 n— 1 

+ - E E ( 5 222/(< / ) n , d 22222 f(Blf n , Bg)) ( e , /B) « fc/B ) fl ( ei/B) tf I/B ) fl 



fe,Z=0 



(4.31) 



(4.32) 
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1 n— 1 

fc,Z=0 



:n) + c (n) + d (n) + e (n) . 



W + b W + c i«; 



Using Lemma I3TT1 (i) and (ii), we have that a (n) , c (n) and d (n) tends to zero as n 
Lemma EH] (iii), we have that tends to zero as n — > oo. Finally, observe that 



Lemma 13.11 (iii 

1 n— 1 

fc,Z=0 



oo. 



Using 



^= £ 5 (w^i, Bg,) ^ 222 /(< , 5$)) ^ /B> ^ + 

- ^ 5 (^2222/(5^, 5^)^222/(5^,5,^)) f (e fc/nJ <5 fe/n ) fl + — ^= J {e l/n , 8 l/n ) 
k,l=0 v v 7 



Therefore, using Lemma l3~Tl (i) and (iv), we have 

2 i 2 

= ^ ^E^22/(5«,5f/j + (1) 
fc=0 



E 



(4.33) 



On the other hand, we have 



E 



r n— 1 1 Ti—l \ 

\ v fc=0 z=o / 

1 n— 1 

E ^(^/(< n ^f;ja 2222 /(5«,5^)A5g)j 
v fc,/=0 
1 n— 1 

-W7fi E S (^222/(5« ,5f/j5 2 222/(5«,5^))( £fe/nj 5 fe/n )^ 
9222/(5^, 5^)522222/(5^, B« )) (s l/n , 8 k/n )„ 



^ n— 1 

v fc,Z=0 
1 n— 1 

= 4^2 E 5(^222/(5^ ,5f;ja 22 22/(5«,5^)) 
fc,Z=0 

n-1 

? fc/n> ^fc/L) -922222/(5^, fijyj) (e, /n , £ fc/n )s 



1 n—1 

v fc,£=0 
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We immediately have that the second (see Lemma 1341 (iv)) and the third (see Lemma I3~T1 (ii)) 
terms in the previous expression tends to zero as n — > oo. That is 



71-1 



71-1 



+ o(l). 



E [\T. ^2/(B« , B$ n )ABg x A d 2222 f(B[)l B\%) 

V V fc=0 (=0 
^ n— 1 

= E 2^^ 92222 f( B k/n> B k/n 
k=0 

We have proved, see (I4~33j) and ([HMD , that 



(4.34) 



71-1 



fc=0 



fc=0 



0. 



This implies (|4~32|) . 

The proof of (I4.2ip follows directly from (14.23P by exchanging the roles played by B^ and 
B^ 2 \ On the other hand, by combining Lemma [3.21 with the following basic identity: 



we see that (|4.27p is also a direct consequence of (|4.32|) . Finally, (I4.26P is obtained from (|4.27|) by 
exchanging the roles played by B^ and B^ 2 \ 



Proof of U-22\ ), U-24\) and M-28\} . By combining Lemma [321 with the identity 



n 



we see that H4.24H is easily obtained through a Riemann sum argument. We can use the same 
arguments in order to prove (I4.22p . Finally, to obtain (|4.28l) . it suffices to combine Lemma l3~2 
with the identity 



Proof o/ (4.29\) and U-SOj) . We only prove (|4.30p , the proof of (|4.29D being obtained from 
(|4.30p by exchanging the roles played by B^ and £?( 2 ). By combining (|4.31l) with Lemma [3.21 it 
suffices to prove that 



1 71 — 1 

-4^ d 1222 f(B^ ,Bf) ) AB« ^ b 0. 

/ / J-"^J V fc/ n ) k/n> k/n k/n 
V k=0 



But this last convergence follows directly from Lemma 1331 Therefore, the proof of (|4.30p is done. 
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Proof of (4-25\) . We combine Proposition 13.31 with the idea developed in the third comment 
that we have addressed just after the statement of Theorem 11.21 Indeed, we have 



k=0 



for = (BW + BW)/y/2 and /? = (B^ - B^)/y/2. Note that (/3,/3) is also a 2D fractional 
Brownian motion of Hurst index 1/4. Hence, using Proposition 13.31 with g(x,y) = g(x,y) = 



fe=0 

S ^ y ^ / 1 S U /(S«,BW)d(^-W0 - + i C d 1X2 ,f{B^\Bf))ds 

n->oo 2 J 4 J 

L = ^ jf 1 dnf{B?\BV>)dW B + ^ jf 1 d^fiBM^ds, 
for (W, W) a 2D standard Brownian motion independent of The proof of (14.25j> is done. 

Proof of the first point (case H > 1/4). The proof can be done by following exactly the 
same strategy than in the step above. The only difference is that, using a version of Lemma [372" 
together with computations similar to that allowing to obtain (|4.32|) . the limits in (|4.21l) - (|4.28l) 
are, here, all equal to zero (for the sake of simplicity, the technical details are left to the reader). 
Therefore, we can deduce (II. 2p by using (I4.20p . n 
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